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ABSTRACT 



Supersonic flow past two oscillating airfoils is 
analyzed in the thesis using an analytical elementary 
theory valid for low frequencies of oscillation. The 
airfoils may have arbitrary stagger angle. This 
approach generalizes Sauer's solution for a single 
airfoil oscillating at small frequencies in an 
unbounded supersonic flow. 

It is shown that this generalization can provide 
an elementary theory for supersonic flow past two 
oscillating airfoils. This aerodynamic tool will 
facilitate the calculation of pressure distribution and 
consequently the calculation of moment coefficient. 
Torsional flutter boundaries are computed. The results 
for the pitch damping coefficient are the same when 
compared with previous analysis. For arbitrary 
frequencies a linearized method of characteristics was 
outlined. 

The elementary theory that has been developed in 
the thesis can be used for flutter evaluation of 
aircraft carrying external stores. The result of the 
thesis is the derivation of the pitch damping 
coefficient which is necessary to predict the flutter 
conditions. 
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I. INTRODUCTION 



Aeroelasticity is the study of the effect of aerodynamic forces on elastic bodies. One 
of the interesting problems in aeroelasticity is the stability of a structure in wind. Since for a 
given configuration of the elastic body, the aerodynamic force increases rapidly with the wind 
speed while the elastic stiffness is independent of the wind, there may exist a critical wind 
speed at which the structure becomes unstable. Such instability could create excessive 
deformation and failure of the structure. A particular problem is the FLUTTER of structures, 
which is a self-excited vibration phenomenon. 

Flutter is caused by the interaction of all three forces: aerodynamic, elastic and inertial. 
Consider as an example a cantilever wing mounted in a wind tunnel (with the root rigidly built 
in). Suppose the wing is deliberately deflected and then released. At sufficiently low wind 
speed, the oscillation that follows this disturbance will quickly die out. However, at some 
higher wind speed one observes a steady oscillation, which maintains itself. This is the critical 
flutter speed for the wing. At higher speeds, the oscillation will be rapidly divergent causing 
structural failure. 

When large external bodies such as engine nacelles, fuel tanks, or electronic warfare 
pods are added to the wing of an aircraft, the dynamic characteristics of an aircraft will be 
changed. Over the last several decades, there has been considerable interest in the calculation 
of the unsteady aerodynamic forces on aircraft. Many of the methods presently in use for 
numerical calculations involve application of a two- dim ensional theory to solve three- 
dimensional flow problems. Most of these methods do not have the ability to model an 
arbitrarily shaped aircraft However, with certain assumptions, it is possible to simplify the 
flow problems to adequately account for the effects of unsteady aerodynamic forces on aircraft 
components. Unfortunately, the transonic flight regime still makes it very difficult to determine 
the unsteady aerodynamic forces because of the nonlinear shock and flow separation effects, 
which occur in this flig ht regime. 

A good example for the requirement of the wing/ stores flutter analysis is the fighter 

aircraft 
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Fighter aircraft are commonly designed for a primary mission, sucH as air superiority, 
which may require few, if any, wing-mounted external stores. However, for increased 
effectiveness and versatility, many secondary missions evolve which are necessary and require 
the use of a variety of external stores. Thus, many combinations of these external stores must 
be carried at various stations on the wings to achieve the complex, multi-role missions required 
by the operational commands. Sometimes, the wing stores are not available to the aircraft 
manufacturer even during the production of an aircraft and in that case, the flutter analysis 
should be conducted using data not available from the original design. 

A typical modem tactical fighter can carry a great variety of wing stores and 
consequendy the total number of possible aircraft/store configurations is huge. Therefore, on 
military aircraft the effect of unsteady aerodynamics is computed for both the clean wings (no 
stores) and the wing with tip missiles. Using these two aerodynamic configurations, the flutter 
analysis is probably done on between 300 and 400 selected wing/ store configurations. The 
cost for the above effort would increase by several orders of magnitude if wing/store 
aerodynamics were also considered for each of these configurations. In many cases, this is not 
necessary since the store aerodynamics has a small effect on die flutter speed. However, there 
are cases where neglecting the store aerodynamics will lead to the overestimation of the flutter 
speed. Therefore, the determination of the cost effectiveness of such analysis is a critical step. 

In this thesis, we first review the most important papers, which have been published in 
recent years to analyze wing-store flutter problems. It is evident from the literature survey that 
most approaches are based on modem CFD (Computational Fluid Dynamics) approaches 
which require extensive computations. 

However, there is still some merit in simplified analyses, which reveals the major 
physical effects. Limiting our approach to low supersonic flight speed we propose to replace 
the multiple stores underneath the wing by a second wing and to apply to a 2-D problem the 
so called “ strip theory approach” which reduces the three dimensional flow problem. This 
simplification has the advantage of accounting for the interference effects between two lifting 
surfaces. Moreover, these interference effects will be stronger than those caused by the actual 
stores. The insights which can be attained with such an analysis, therefore will represent an 
extreme limiting case, but it will nevertheless reveal important physical effects. 
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The importance of such interference effects is well known in transonic and supersonic 
wind tunnel testing and in the operation of transonic and supersonic compressors. Platzer 
(1973) presented an analytical theory for the analysis of oscillatory supersonic wind tunnel and 
blade interference effects. This theory is based on the assumption that the actual airfoil or 
compressor blade can be replaced by a flat plate and that it oscillates with low frequency. This 
makes it possible to expand the solution in powers of the reduced frequency and to retain only 
the zeroth and first order term. Platzer (1973) showed that this theory provides a convenient 
analytical way to estimate the pitch damping as a function of supersonic Mach number and 
pitch-axis (elastic axis) position. 

In this thesis, we modify Platzeris solution to the case of two airfoils, which are in 
close proximity to each other. Two cases need to be distinguished, namely the one with 
supersonic leading edge locus and the case with subsonic leading edge locus. These two cases 
and the basic theory will be explained in chapter III. The extension of Platzer’s theory will be 
given in chapter IV. This is followed by a chapter on the method of characteristic approach, 
which can be used to compare with the analytical results. The final chapter is devoted to the 
flutter analysis. 
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II. LITERATURE SURVEY 



EVALUATION OF METHODS FOR PREDICTION AND PREVENTION OF 
WING/STORE FLUTTER. 

Authors): Pollock, S. J. Sotomayer, W. A.; Huttsell, L. J.; Cooley, D. E 

Source: Collect Tech Pap AIAA ASME ASCE AHS Struct Dyn. Mater. Conf. 22nd AIAA 
Dyn. Spec. Conf., April 6-10 1981, Atlanta, GA, New York, NY, pp. 362-372 

In response to the need to reduce costs and improve safety for flutter evaluation of 
aircraft carrying external stores, the Flight Dynamics Laboratory (FDL) has sponsored several 
efforts in the technical areas of unsteady aerodynamics, flutter prediction, and active flutter 
suppression. This paper discusses each of these three areas as they relate to wing/store flutter 
and presents specific examples from analyses and tests. Steady and unsteady pressure 
measurements were obtained in a wind tunnel at subsonic, transonic and supersonic speeds on 
a fighter wing, tip-mounted launcher and store, and underwing pylon and store. Store flutter 
calculations were performed using both calculated and measured data to determine the 
influence of store aerodynamics on the flutter characteristics. To improve the accuracy and 
reduce the time and costs of flutter evaluations on the many store configurations carried by Air 
Force fighters, the Flight Dynamics Laboratory (FDL) has sponsored several programs in the 
technical areas of unsteady aerodynamics and flutter predictions. Also, the FDL has been 
sponsoring several efforts to explore the potential of active flutter suppression systems using 
feedback control techniques to provide the required stability and to avoid speed placards. 

This paper reports on some of the FDL research related to wing/store flutter 
prediction and prevention. The research includes an unsteady aerodynamic measurement 
program for a representative fighter wing, with and without tip missile and underwing store 
with test data covering the Mach number range 0.6 to 1.35. Flutter analyses based on an FDL 
computer program specifically for use on aircraft with external stores are also described. 
Flutter trends using this computer program are presented for the wing with and without stores 
based on sectional force coefficients from wind tunnel measurements and from theoretical 
calculations. A brief description is given of some FDL programs in active flutter suppression, 
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and typical results are presented for wings with stores, which indicate significant potential for 
improvement in flutter speeds. 

Flutter trends for three different fighter wing/ store configurations were calculated 
using the FACES flutter analysis procedure. Use was made of available modal vibration data 
and measured aerodynamic data. The predicted Mach number trends gave minimum flutter 
speeds at transonic and low supersonic speeds as would be expected from the trend of 
measured center of pressure and lift-curve slope data. Although the aerodynamics on the tip 
launcher had a somewhat detrimental effect on flutter, the aerodynamics on the tip store had a 
much larger detrimental effect. For the underwing store, the aerodynamic effects were 
beneficial for flutter since the pylon had the effect of decreasing the aerodynamic loading on 
the outer portion of the wing. 
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WIND TUNNEL TESTS ON A FIGHTER AIRCRAFT WING/STORE FLUTTER 

SUPPRESSION SYSTEM 



Authors): Turner, M. R 

Source: AGARD, Rep. on a Coop. Program on Active Flutter Suppression, Paper Presented at 
the Structural and Material Panel Meet., 50th April 1980, Athens, Greece AGARD. 

A special system designed using analytical data was tested on a YF- 17 aircraft model in 
the NASA Langley 16' wind tunnel and succeeded in meeting the requirement to increase the 
flutter dynamic pressure by 70% at Mach equals 0. 8. The system was designed using a novel 
procedure, which provides these stability margins, uses minimum control surface movement in 
turbulence, and can be designed using either analytical or empirical data. Two wing tip 
accelerometers and a leading edge control surface were used. Empirical open loop transfer 
during the test showed that the analytical data overestimated the response 
to leading edge control surface excitation. 



functions obtained 
of the flutter mode 
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AERODYNAMICS USING TIME-LINEARIZATION TRANSONIC FLUTTER 

ANALYSIS 



Authors): Wong, Y.S., Lee, B.H.K.; Murty, H.S. 

Source: Journal of Aircraft, v. 30, January-February 1993, pp. 144-145 

A survey of the progress made in the development of numerical simulation techniques 
for unsteady transonic flow calculations are presented. Computational methods in three- 
dimensional unsteady transonic flows concentrate mainly on the transonic small disturbance 
equation and time-lineraization approach. An algorithm is introduced for solving flutter 
occurrence. 
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WING FLUTTER BOUNDARY PREDICTION USING UNSTEADY EULER 

AERODYNAMIC METHOD 



Authors): Lee-Rausch, Elizabeth M.; Batina, John T. 

Source: Journal of Aircraft, v. 32, 2 March- April 1995, ALAA, Washington, DC, USA, pp. 416- 
422 



Modifications to an existing three-dimensional, implicit, Euler/Reynolds-averaged 
Navier-Stokes code for the aeroelastic analysis of wings are described. These modifications 
include the incorporation of a defor min g mesh algorithm and the addition of the structural 
equations of motion for their simultaneous time-integration with the governing flow equations. 
This article gives a brief description of these modifications and presents unsteady calculations 
that check the modifications to the code. Euler flutter results for an isolated 45-deg swept- 
back wing are compared with experimental data for seven freestream Mach numbers that 
define the flutter boundary over a range of Mach number from 0.499 to 1.14. These 
comparisons show good agreement in flutter characteristics for freestream Mach numbers 
below unity. For freestream Mach numbers above unity, the computed aeroelastic results 
predict a premature rise in the flutter boundary as compared with the experimental boundary. 
Steady and unsteady contours of surface Mach number and pressure are included to illustrate 
the basic flow characteristics of the time-marching flutter calculations and to aid in identifying 
possible causes for the increase in the computational flutter boundary. 



9 



APPLICATION OF TRANSONIC SMALL DISTURBANCE THEORY TO THE 

ACTIVE FLEXIBLE WING MODEL 



Authors): Silva, Walter A.; Bennett, Robert M 

Source: Journal of Aircraft, v. 32, January-February 1995, AIAA, Washington, DC, USA, pp. 
16-22 



A code, developed at the NASA Langley Research Center, is applied to the active 
flexible wing wind-tunnel model for prediction of transonic aeroelastic behavior. A semispan 
computational model is used for evaluation of symmetric motions, and a full-span model is 
used for evaluation of antisymmetric motions. Static aeroelastic solutions using the 
computational aeroelastidty program-transonic small disturbance are computed. Dynamic 
(flutter) analyses are then performed as perturbations about the static aeroelastic deformations 
and presented as flutter boundaries in terms of Mach number and dynamic pressure. Flutter 
boundaries that take into account modal refinements, vorddty and entropy corrections, 
antisymmetric motions, and sensitivity to the modeling of the wingtip ballast stores are also 
presented and compared with experimental flutter results. 
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WING-STORE FLUTTER ANALYSIS OF AN AIRFOIL IN INCOMPRESSIBLE 

FLOW 



Authors): Yang, Zhi-Chun Zhao, Ling-Cheng 

Source: Journal of Aircraft, v. 26, 6 June 1989, pp. 583-587 

The flutter of two-dimensional airfoil with external store is analyzed to investigate the 
effects of pylon stiffness on flutter speed Among the 40 configurations studied, five were 
tested in the wind tunnel to verify the analytical results. The variations of wing-store flutter 
speed with the pylon stiffness can be divided into three types. The curves of the normal and 
flutter frequencies vs pylon stiffness have the same pattern. They can be sketched 
approximately by the aid of the normal frequencies of the two degenerated two-degree-of- 
freedom systems, i.e., and the freely hinged and rigidly connected store cases. A limiting flutter 
speed for very small pylon stiffness is deduced, which is useful to identify which type of flutter 
the configuration studied belongs to. 
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STUDY OF THE EFFECT OF STORE AERODYNAMICS ON WING/STORE 

FLUTTER. 



Authors): Turner, C. D. 

Source: Collect Tech. Pap. ALAA ASME ASCE AHS Struct Struct Dyn. Mater. Conf. 22nd 
AIAA Dyn. Spec. Conf., April 6-10, 1981, and 1981, Atlanta, GA, AIAA (CP 812), New York, 
NY, pp. 343-351, Paper: 81-0604 

This study represents the first systematic analytical study of the effect of store 
aerodynamics on wing/ store flutter. A large number of wing/store single carriage 
configurations and parameters were included in the study; multivariate analysis techniques 
were used for the first time to analyze wing/store configurations, modal data, and flutter 
results. The results of the multivariate analysis indicate that it may not be possible to develop 
general guidelines, but it is possible to develop specific guidelines for use with a particular 
aircraft This study was the first attempt to do a systematic analytical study of the effect of 
store aerodynamics on wing/store flutter. To determine this effect flutter analyses were done 
on four aircraft with single carriage of three basic store types. In all 308 configurations were 
analyzed with and without store aerodynamics (tip missiles, tip tanks, and underwing stores). 
The effects of stores have been analyzed. The results of the factor analysis indicate that it may 
not be possible to develop general guidelines, but it is possible to develop specific guidelines 
for use with a particular aircraft The conclusions of the study, as far as it concerns the effect 
of store aerodynamics on wing-store flutter, show that 60% of the tested configurations 
require a change in the flutter speed no more than 7%. Therefore, there is no need to evaluate 
the requirements for aerodynamic modeling of the store. However, for the rest of the 
configurations there is such a need. Moreover, 75% of the last reveal that nonrealistic flutter 
results are obtained without taking into consideration the store aerodynamic. Of course, the 
results of this study represent the configurations that were used to generate the data base. 
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AERODYNAMIC MODELING OF AN OSCILLATING WING WITH EXTERNAL 

STORES 

Authors ):, Sotomayer W. A., Dusto A. R, Epton MA., Johnson F. T. 

Source: AIAA , New York, NY, pp. 243-252, Paper: 81-0609 



An analysis of the steady and unsteady aerodynamic forces acting on a fighter aircraft 
wing with stores has been done. Computations were performed with paneling methods 
capable of presenting arbitrary aircraft configurations in subsonic and supersonic flow. 
Interference effects from a tip store and an underwing pylon store in varying stages of 
completion were also analyzed Detailed comparisons between experimental data and 
numerical computations are also made. 

In order to calculate the pressures, forces and moments on an aircraft a numerical 
approach is proposed. Numerical computations are based on solving an integral equation 
formulation of the flow problem being considered. Steady state calculations were done with 
two numerical methods. The first of these is a pilot code developed by Johns. In this method, 
distributions of linearly varying sources and quadratically varying doublets are used to 
represent the aerodynamic surfaces. Wakes are represented by doublets with constant 
streamwise strength and linear variation in the spanwise direction. A method developed by 
Woodward was also employed in making numerical calculations. This method utilizes 
distributions of sources and vortices to represent an aerodynamic surface and the wakes shed 
from various components. 

For unsteady flow, each of the methods makes use of sources and doublets as the basic 
aerodynamic singularities. In a method developed by Johnson for unsteady subsonic flow, 
distributions of sources and doublets are used to represent the aerodynamic surfaces. In the 
wake, doublets of fixed strength and position are used to represent the effects of unsteady 
wake motion. Calculations were also done with a special doublet lattice method. In this 
method, a wing is represented as a sheet of doublets, components such as the fuselage or a 
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store or nacelle can be represented with an axisymmetric distribution of sources with doublets 
on the surface. 

Approximations are imposed on these distributions and on the kernel functions used to 
represent their influences on the flow. These approximations are briefly described along with 
the flow boundary conditions imposed on the aerodynamic surfaces and wakes. 

It is interesting that the small disturbance partial differential equation is being used for 
unsteady in viscid incompressible flow along with the linearized form expressed in terms of the 
velocity potential cp. In addition, simple harmonic motion is assumed as in the thesis. 

<p=Re[<p*e^] where <p* is a complex quantity. 

Moreover using the Helmoltz equation a more simple relationship is established for (p* with 
the form of an integral equation. An approximate solution of the last equation permits an 
evaluation of the entire flowfield. Numerical solutions of the last equation are achieved 
through an aerodynamic influence coefficient method. The boundary conditions are 
established taking into consideration the configuration of the F-5 aircraft, a relatively old 
aircraft. The aircraft is comprised of the following components: 

1) wing 

2) underwing pylon 

3) missile rack 

4) missile body. 

As in the thesis, the above surfaces are approximated as infinitesimally thin surfaces so 
the thin wing theory applies except the for missile body. The boundary conditions are 
expressed in terms of Taylor series expansion about reference surfaces taking into 
consideration the thickness. 

Initially boundary conditions for steady state flow are made use of and unsteady 
boundary conditions are taken as a perturbation about a steady state condition. In addition, 
variations in Mach number, frequency of oscillations, and interference effects arising from 
component build up are presented. 
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For the wing missile, the boundary conditions are of mixed type (Dirichlet and 
Newman). On the surface of the missile, the boundary conditions are similar to the other parts 
of the aircraft. 

Using a Woodward method, a lifting surface is divided into numbers of aerodynamic 
panels number, containing distributions of sources and vortices. A source distribution 
represents a fuselage [or a pod or external store]. Wing thickness is represented by a linearly 
varying source distribution in which the strength is equated to the chordwise slope of the wing 
thickness. Camber, twist, and lifting effects are represented by a linearly varying vortex 
distribution where the strength is determined to satisfy tangential flow at panel control points. 
An iterative procedure is employed in solving the boundary value problem. For analysis 
problems, the surface slope is described, and singularity strengths are determined by inverting 
the matrix of aerodynamic influence coefficients. With the strengths of the aerodynamic 
singularities known, the u, v, and w velocity components at a given point may then be 
determined. Pressures, forces, and moments are calculated by numerical integration. 

In the doublet lattice method, the aerodynamic surface is subdivided into a series of 
infinitesimally thin panels. Along the quarter chord line of each of these panels is contained a. 
distribution of acceleration potential doublets. The strength of the distribution is constant but 
is not known. Specification of the normal velocity at a set of points on the surface determines 
the loading on each element Locating the lifting elements at the quarter chord and the 
collocation point of the three-quarter chord of the midspan of each panel usually results in 
reasonable success. An approximation for the lifting pressure coefficients in terms of the 
induced normal velocities may then be found. lift, moment, roll and generalized forces may 
then be calculated. 

Recently, a series of wind tunnel tests were conducted by NLR of the Netherlands with 
sponsorship from the Air Force. These tests involved measurements of steady and unsteady 
pressures and forces on a model of the F-5 wing with external stores. Several external store 
arrangements, which are described and shown in this paper, are compared with numerical 
calculations. Experimental results and numerical calculations show that 

1) For subsonic flow, the interference due to the tip store is much greater on the 
outer portion of the wing than it is on the inner portions. 
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2) The interference effect due to the pylon store is significant on the lower surface of 
the wing for subsonic flow; there are also noticeable effects on the upper surface of the wing at 
subsonic speeds. 

3) For supersonic speeds, the zone of influence remains localized within the Mach 

cone. 

In general, numerical calculations and experimental results agreed reasonably well with 
each other. 
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III. PROBLEM FORMULATION 



The interference between the airfoil (wing airfoil) and the wing stores is examined by 
making certain assumptions that simplify the problem. The wing airfoil is replaced by an 
oscillating flat plate airfoil and the wing stores by another oscillating flat plate in distance d 
from the wing airfoil Two cases need to be distinguished, namely; 




Fig. 1 

1) Supersonic flow past two oscillating airfoils with supersonic leading edge locus 
shown in Fig. 1 (above) 

2) Supersonic flow past two oscillating airfoils with subsonic leading edge locus 
shown in Fig. 2 
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The differences are obvious and will be further analyzed during the examination of 
each case. In both cases multiple reflections occur. As explained in the preface the interference 
between the two airfoils under certain conditions could result in unsteady flow phenomena like 
flutter. 

Therefore, the study of the oscillating airfoils is important in predicting the flutter 
characteristics or dynamic response of the wing airfoil when different stores are attached. For 
the case shown in Fig. 1: 

tan/? < cota (Ill-la) 

Where the case shown in Fig. 2 requires: 

tan/? > cota (M-lb) 

The flow is assumed a non- viscous compressible two-dimensional flow of a perfect gas 
governed by the continuity equation, the Euler, and the energy equation. The continuity 
equation is: 



Dp du dv 

— — + p — + — 
Dt dy j 



= 0 



The Euler Equations: 



El + L& = 0 

Dt p dx 



Dt p dy 



The Energy Equation: 




Dt 



(III-2) 



(III-3) 



(IH-4) 



(HI-5) 
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Because of the assumption of small amplitude oscillations of the airfoils, all flow quantities are 
considered small perturbations linearly superimposed on free stream quantities. The velocities 
can therefore be written as: 



U — ll x + u 


(IH-6) 


c = c +c 

oo 


(HI-7) 


V = V 


(in-8) 



In addition, the pressure and density perturbations are linearly superimposed on freestream 
quantities. Therefore: 



£ 

II 

1 

8 


a n-9) 


Ap = p-p x 

The local velocity of sound is given by: 


(III-10) 


- f dp'? 

° l dp) S=Ct 

Furthermore: 


(III-ll) 


p 

— = const 
P r 

Taking the total differential: 


(III-12) 


s? \|- 
1 

•fc - 

V 

II 

o 


(in- 13 ) 
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(ni-14) 




So: 




(111-15) 



From AA 4318 class notes, the surface boundary conditions will be calculated Requirement is 
that no flow will penetrate the solid surface of the airfoil resulting in the flow tangency 
condition: 

(y~V s )-n - 0 (UI-16) 

where V,V S are the velocities of the free stream and the surface over which the fluid flows.If 
the body is described by the equation: 

F(x,y,t) = 0 (III-17) 



then the normal to the surface at any point is given by the following relation: 



>F| 



(in-i8) 



Taking into consideration that the substantial derivative of a surface particle should be zero: 
dF - 

— + V S ‘VF = 0 (III-19) 

dt 

one obtains from (11-16): 

FVF = V s • VF (III-20) 

and therefore: 
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(III-21) 



DF 

Dt 



dF — 

= —— + V -VF = 0 

dt 



Or for two dimensional flow: 



DF dF dF dF 

= h it — + v — = 

Dt dt dx dy 



(Til-22) 



Particles in contact with the surface must have the same normal velocities as the surface. This 
means that the rate of change of F is zero hence: 



DF 

Dt 



= 0 



(III-23) 



For an airfoil the equation of the upper surface can be written as follow: 

F u (x,y,t) = y-y u (x,t) = 0 (III-24) 



Where y u is the distance from the chord line to the upper surface. The equation of the lower 
surface can be written in a similar way: 



F L (x,y,t) = y-y L (x,t) = 0 



(III-25) 



Aty=y u : 



DF u _ d}’u ( -x d}’u 

Dt dt K dx 



+ v = 



0 



Aty=y L : 



DFt dy L - dyi 

— — = — — - (it) ■ -F-F- + v = 0 

Dt dt dx 



(III-26) 



(III-27) 



tyu _ &L _ { 
dy dy 



(III-28) 
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Thus, the normal flow velocity can be written: 



v = + (III-29) 

dt dx 

Applying the assumption of linear perturbation theory one has u = u x +u , therefore: 

v = ^+ u „^ + u ^ 7 ..y = y u (x,t) (HI-30) 

dt ox ox 



V = Ql + u x ^- + u^,..y = y L (x,t) 
dt “ dx dx y yLX 1 



(111-31) 



Neglecting the higher order terms one obtains: 



v = %-+ ^->y = y u (x,t) 

dt dx 



(III-32) 



dt * iV ' 



(III-33) 



This normal velocity will be expanded as Taylor series of the normal flow velocity around 
y=0: 

v(,, yi> ,) = v(,,0V) + * + (.11-34, 

oy 21 a y 

v ( x > y L » 0 s = v ( x » 0+ > r )+Ti ^ a h° + ail-35) 

qy 2! a _y 

Using the assumption of thin airfoil and small linear perturbation theory the higher order 
terms may be neglected and one obtains: 



v (*>T u ,')= v (*>°V) 



(m-36) 
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(m-37) 



Therefore: 



v[xfi\,)=^ +u ^,..y = 0- (01-38) 

v[xfi-,t)=^+u.^,y = 0- ' (01-39) 

at ox 

For simple harmonic motion of the airfoil: 

y. - KW“ (io-40) 

y L = (10-41) 

Therefore, the Unearned flow tangency equations can be written: 

v(x,0V)= [icoh u {x) + u„^]e“*,..y = 0 + (HI-42) 

3x 

v(x,0v)=[/a !Vl ( J r) + I /„^)e",..y = 0- ' (10-43) 



The above relationships describe the normal velocity on the lower airfoil (store airfoil). 
Similarly, one obtains for the upper airfoil (wing airfoil) at y=d: 

v(x,d~,t)= [icoh L {x) + uj^e^.y = d~ (III-44) 

ox 

if the store airfoil executes a pitch oscillation about x = x 0 one has: 
y = -e(x-x o y rt (HI-45) 

v(jr,0 + ,0 = [-uj - icoO{x- x^)[e iaM at y = 0 + (III-46) 
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and similarly for the wing airfoil: 



y = -0(x-B-x o )e i(wt+S) (III-47) 

v(x,d~,t ) = [- 0{u x cos 8- 6>.sin£(jt - B - ;c 0 }- i6{u x sin# + cy.cos#(;t - B - * 0 )}]e ,a " 

(III-48) 
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IV. ELEMENTARY THEORY 



The case of supersonic leading edge locus will be examined assuming first that both 
airfoils are oscillating slowly in a supersonic flow. The airfoils execute vibrations of identical 
modes and amplitude but with different phase angle 8 between them. It is obvious that the 
flow field is such that no disturbance can propagate upstream of the airfoil leading edges. 
Consequently, the perturbation potential can be written as: 

<I>(x,y,t) =(p(x,yh)e to (TV-1) 

where x, y are non dimensional coordinates and t is non dimensional time. The differential 
equation for the perturbation potential in a non-dimensional form then is: 

cota-cp :a -(p yy +2ikM 2 (p x -k 2 M 2 (p=0 (IV-2) 

The linearized boundary conditions are: 

<p y (x,0)=v(x) on the lower(store airfoil) at y=0 (TV'S) 

where v(x) is the downwash amplitude. 

'5 

(p y (x,d)=— e 1 v(x— dtanP) on the upper(wing) airfoil at y=d (IV -4) 

where it is assumed that the two airfoils oscillate with a phase angle 8. 

Assuming oscillation around Xq v(x) becomes: 

v(x)=0[l +ik(x— X q)] and without loss of generality we can set 0=1 (IV-5) 

For small oscillations (low frequency) the potential cp can be expanded and neglecting the 
higher order terms of k one can write: 

<p(x,y,k) =X(x,y) +kT(x,y) (IV-6) 
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Therefore the differential equation (II-2) becomes: 



cot 2 a*X D -X yy =0 (TV-7) 

cota- v F yy -T yy + 2iMX=0 (IV-8) 

A general solution of the above system of equations has been given by Sauer (1950). It reads: 

X(x,y)=g(z) (IV-9) 

'P(x,y)=h(z)+iMyg(z)/cosa where z=x-ycota (TV-10) 

v F(x,y)= h(z) +iMy g(z)/cosa where z=x+ ycota (TV-11) 

h(z) and g(z) are arbitrary’ functions for positive arguments of z and zero for z<0, it is obvious 
that the expressions of the above functions should be consistent with the last assumption. 

Using the two solutions for left and right running Mach waves in the supersonic flow 
field of slowly oscillating airfoils pressure coefficients can be obtained. In order to do so the 
flow field between the two airfoils has to be divided into several zones. The number of zones 
depends upon A. It is obvious that for A>1 where A=dcota there is one zone along each 
airfoil, for 0.5<A<1 there are two, for 0.33<A<0.5 there are three and so on. In the first zone 
I there is no interference from the airfoil wing (upper wing), which means that the lower airfoil 
(zone I only) does not sense the upper airfoil In that case, the solution has the form: 

9 ( x »y>k) =g(z) + k { h (z)— iMyg(z)/cosa} 

The boundary conditions that must be satisfied at y=0 are: 

X y =-cota-g'(z)=— 1 

T y = -cotah' (z)-iMyg(z)/ cosa=-i(x -Xq ) 

From the last two equations it is obvious that: 
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g(z)=ztanct 



h (z) =— iztana(x 0 + ztan 2 a/2) 

Since we are assuming a slowly oscillating airfoil, the solution to the PDE has the form: 
cp (x,y,k) =g(z) + k [h (z)-iMyg(z)/ co sa] 
where z=x— ycota 

In order to calculate the pressure distribution on the upper surface of the lower airfoil in the 
zone I the following relation should be used: 

cp=-2fikcp +cp i ] 

Substituting the known functions cp and cp, 

cp=ztana+k{— iztana(x 0 +ztan a/2) — (iMyztana/cosa)} 

cp x =d/dx[ztana +k{-iztana(x 0 +ztan"a/2) — (iMyztana/ cosa) } ] = 

= tana— iktana [xq+ ztan“a- (My z/ cosa)] = 

=tana-ik[tana(x 0 +ztan 2 a)+( y/ sinacosa)(sina/cosa)]= 

2 2 

=tana-ik[tana(x 0 +ztan“a)+fy/cos‘'a)] = 

• 2 2 
=tana-ik[tana(x 0 +(x— ycota)tan a+(y/cos a)]= 

=tana-ik[tana(x 0 +x tan a— ytana+(y/cos a)]= 

2 

=tana-ik[tana(x 0 + x tan a— ycota] = 

Therefore, the pressure distribution is given by the following relationship: 
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q)=-2{tana-ik[tana(x 0 +xtaa 2 a-ycota]+ik[ztana+k{-i2tana(x 0 +ztan 2 a/2)-iMyztana/ 

cosa}]} 

Taking into consideration that higher order terms of k are neglected the pressure distribution is 
given by the following relationship: 

cp=2tana{— 1+ikjxg+x (tan 2 a— l)+2ycota]} 

and therefore on the airfoil fory=0: 

cp=2tana{-l+ik[xo+x (tan 2 a-l)]} 

In the zone I* the perturbation potential can also be written directly from Sauer* s single airfoil 
solution: 

<P( x >y>k) =gi(z 1 )+k[hi(z 1 )+ iMyg(zi)/cosa] 

Where: 

Z!=x-A-B +ycota 

The boundary conditions that must be satisfied at y=d are: 

X y =cota-g'(z 1 ) = -e 5 

% =cota-h 1 '(z 1 )-iMg 1 (zi)/cosa +iMy g/(z 1 )=-ie , ^(x-x 0 -B): 

After integration and algebraic manipulations: 

gi( z i)--Zitana-e^ 

hi(zi)=iz 1 tana-e‘°(x 0 +z 1 tan 2 a/2+Md/cosa) 
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where again the condition has been imposed that the functions gi(zj) and hi(zi) are zero for 
negative and zero arguments. 

Following the same procedure as in zone I, the perturbation potential (p will be 
calculated substituting the known functions gi(zi) and h^z^: 

S ’8 "S 

(p=-Zitanae 1 +ik{z 1 tana e‘°(x 0 +z 1 tan 2 a/2+Md/cosa)— Myzjtanoce /cosa} 



"8 *8 

<p x =— tanae 1 +iktana e'° { (xo+z^a^a+Md/cosa)— My/ cosa} = 



’8 *8 0 
r — tana-e 1 +iktana-e ,0 {(xo+(x-A— B+ycota)tan“a+Md/cosa)-My/ cosa} ; 



'8 '8 o 

=-tana-e' +iktana-e ,0 {(xo+(x— B)tan 2 a-A tan 2 a+ycota+d/sina-cosa) 



-y/sinacosa) = 



8 *8 

= — tana-e l +iktana-e'° { (xq+(x— B)tan 2 a+ A— ycota} 



Because d=A/ cota , tana-1 /cosa-sina=cota 

This gives the following result for the pressure distributions on the lower surface of the upper 
airfoil (wing airfoil) in zone I*: 



cp— 2tana-e'^— 2iktana-e 1 ^ { (xq+ (x— B) (tan 2 a— 1) +2A— 2ycota} = 

.C 

cp=2tana*e' { 1— ik{ (xo+(x-B)(tan 2 a-l)+2A— 2ycota} 
and therefore at y=d: 

cp=2tana-e'^ { 1 — ik { (xo+ (x-B) (tan 2 a— 1) } 
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Assuming that 0.5<A<1 in addition to zones I, I* there are two more zones II, II*. In 
zone II the perturbation potential is due to the initial waves from zones I, I* plus the reflected 
wave at zone II.So the reflected waves must be determined by adding the functions g>(z£ and 
h 2 (z 2 ). In the zone II Sauer’s solution for a single airfoil applies and following the same 
procedure as previously: 

(p(x,yh) =g(z)+k[h(z)-iM(yg(z)/cosa l+gtCzO+kfh^zO+iMfygifo)/ cosa]+ 
g 2 (z 2 )+k(h 2 (z 2 )-iMyg 2 (z 2 )/cosa] 

The variable z has the following form: 
z 2 =x-2A+ycota 

This ensures that z 2 <0 upstream of the leading edge of the airfoil. 

The boundary conditions that must be satisfied at y=d are: 

.C 

X y = -cota-g' (z) +cota-g 1 '(z 1 ) +cota*g2' (z 2 )=-e l 

r 

v J , y = -cota-h 1 '(z)+cota-h 1 '(z 1 )+cota-h 2 '(z 2 )-iM[g(z)-gi(z 1 )-g 2 (z 2 )]/cosa+iMdcota[g 1 (zj) 

ft 

+ gi ( z i)+g2 (z 2 )]/cosa=-ie ,0 (x-xo-B) 

After integration and algebraic manipulations s imil ar to the previous calculations for zones I, 
I*: 

g 2 ( z 2 )=ztana 

h 2 (z 2 ) =— iz 2 tana (xq+ z-J 2tan 2 a +2Md/cosa) 

Therefore, the perturbation potential is given by the following relationship: 
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.g ,g 

<j> I =2tana-tana-e‘ -iktana(xo+z 2 tan 2 a+2Md/ cosa)+iktana-e‘°(x 0 +Zitan 2 a+Md/ cosa) 

2 "S 

— iktana(xo+ ztan a )-ikMy [tana-tana + tana-e' ]/cosa= 

=2tana— tana-e' — iktana[(xo+z 2 tan 2 a+2Md/cosa)— e (x 0 + Zjtan 2 a +Md/ cosa)— ( xq+ ztan 2 a 
)+Mye'^ /cosa]= 

,g o 

=2tana— tana-e* — iktana[xo+(x— 2A+ycota)tan 2 a+2Md/cosa-e (x 0 +(x-A-B+ycota)tan 2 a 
+Md/ cosa)-(x 0 + ( x— ycota) tan 2 a)+Mye‘°/ cosa] = 

,g g ^ 

=2tana— tana-e' -iktana[x 0 +(x-2A)tan“a+2Atana/sina-cosa-e' (x 0 +(x— A-B+ycota)tan‘'a 
+Atana/ sina-cosa)— ( xq+x tan 2 a-ye'°/sina-cosa)]= 

=2tana— tana-e* -iktana[x 0 +xtan 2 a+2A-e‘°(x 0 + (x-B+ycota)tan 2 a+A)+x 0 +xtan 2 a+ye* 
/sina-cosa] = 

=2tana— tana-e' — iktana[2x 0 +2xtan“a+2A-e‘°(x 0 +(x— B+ycota)tan 2 a+A)+ye /sina-cosa]= 
=2tana-tana-e*^— iktana[(2-e'^ )(x 0 +xtan“a+A)+ e‘^(Btan"a+ ycota)] = 

=tana{2— e'^— ik[(2— e'^ )(x 0 +xtan“a+A)+e'^(Btan”a+ ycota)]} 

This gives the following result for the pressure distributions on the lower surface of the 
upper airfoil (wing airfoil) in zone II: 

Cp=-tana{2-e‘^+ik[e'^(x 0 +(x-B)(tan 2 a-l)+2A-2ycota)- 
2(x 0 +x(tan 2 a- 1)+2A)} 
and for y=d we have: 
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Cp — tana{2-e‘^+ik[e l ^(x 0 +(x-B)(tan 2 a— l))-2(x 0 +x(tan 2 a-l)+2A)} 



In zone IP the perturbation potential is due to the initial waves from zones I, I* plus 
the reflected wave at zone II, IP.So the reflected waves must be determined by adding the 
functions gj(z 3 ) and h 3 (z 3 ). Therefore, in the zone IP Sauer’s solution for a single airfoil 
applies and following the same procedure as previously. 

<P(x»y>k) =g(z)+k[h(z)-iM(yg(z)/ cosa J+g^zO+kfh^z^+iM^g^Z!)/ cosa] + 



g 3 (z 3 )+k[h 3 (z 3 )— iMyg 3 (z 3 )/ cosa] 

The variable z 3 has the following form: 
z 3 =x-A-B-ycota 

The boundary conditions that must be satisfied at y=0 are: 

Xy =- cota*g'( z) +cota-g/( z x ) -cota-g^ z 3 )= -1 

'By=-cota-h'(z)+cota-h 1 '(z 1 )-cota-h 3 , (z 3 )-iM[g(z)-g 1 (z 1 )+g 3 (z 3 )]/cosa=-i(x-xo) 

After integration and algebraic manip ulations similar to the previous calculations for zones I, 
P: 

g 3 ( z 3 )=— Zjtana-e* 5 

h 3 (z 3 ) =iz 3 tana (xo+z 3 tan 2 a/2+Md/cosa)e i5 

Therefore, the perturbation potential is given by the following relationship: 

5 o ,5 

(p x =tana-2tana-e l0 +iktana{e ld (xo+z 1 tan 2 a+Md/ cosa)+e‘ (x 0 +z 3 tan 2 a+Md/ cosa)— ( xq+ 
ztan 2 a )— My/cosa} 
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,o o o 

cp 1 =tana-2tana-e 1 +iktana{e'°[xo+(x— A-B+ycota)tan‘a+Md/cosa]e [xo+ 
(x-A-B-y co ta) tan 2 a +Md/cosa] — (xo+ (x-ycota)tan 2 a )— My/cosa} = 



<p I =tana-2tana-e 1 +ik{e lO 2tana[x 0 +(x-A-B)tan*'a+Md/cosa]-tana(x 0 +(x— ycota)tan“a)- 
y/cos 2 a} 

(p I =tana-2tana-e 1 +ik{e l0 2tana[xo+(x— B)tan 2 a-Atan‘'a+Adtana/sina-cosa]— tana(x 0 + 
(x-ycota)tan 2 a )-y/cos 2 a} = 

(p I =tana— 2tana-e‘ +ik{e l0 2tana[xo+(x— B)tan“a+A]— tana(x 0 +(x— ycota)tan~a)— 
y/cos 2 a} = 

(p I =tana-2tana-e , °+ik{e l0 2tana[xo+(x-B)tan‘a+A]-tana(x 0 +xtan 2 a+ ycota)} 

Because: y/ cos 2 a+ ycotatan 3 a=— y 

In order to calculate the pressure distribution on the upper surface of the lower airfoil in the 
zone II* the following relation should be used: 

cp =-2fikcp +cp x ] 

Substituting the known functions cp and <p x 

=2tana-4tana- e— 2ik { eztana [xo+ (x-B) tan 2 a + A] - tana (xq+ xtan 2 a + ycota) } 

-2ik { ztana-Zjtana- e — z 3 tana* e'°}+higher order of k terms which will be neglected 
Therefore: 

Cp=2tana-4tana- e 1 ^— 2ik { e*^2tana [xq+ (x-B)tan 2 a+ A]— tana(x 0 +xtan 2 a+ yco ta) } 
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‘8 iS 

-2ik{(x-ycota)tana-(x-A-B+ycota)tana- e— ( x— A-B— ycota)tana- e' } 

After several relatively simple algebraic manipulations: 

Cp=— 2tana[ 1 — 2e'^-ik { 2e 1 ^ [xq+ (x— B)(tan 2 a-1)+2A]— (x 0 +x(tan 2 a-l)+2ycota)}] 

Since the pressures have been calculated for the flow field between the airfoils a similar 
approach will be utilized to obtain the pressures on the bottom of the lower airfoil. In that 
case, the solution has the form: 

(p(x,y,k) = g(z) +k{h(z) +iMyg(z)/cosa} 

The boundary conditions that must be satisfied at y=0 are: 

Xy =cota-g'(z)= -1 

T y = cota-h'(z)+iMg(z)/ cosa=-i(x —Xg ) at y=0 
From the last two equations it is obvious that: 
g(z)=ztana 

h(z)=-iztana(x+ztan 2 a/2) 

The solution again has the form: 
cp(x,yjs)=g(z)+klh(z)+ iMyg(z)/cosa ] 
where z = x+ycota 

In order to calculate the pressure distrib ution on the upper surface of the lower airfoil in the 
zone I the following relation should be used: 

cp =-2[ik(p +cp x ] 
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cp — -z tana + k { — i z tana(x + ztan'a/2) + (iMy ztana / cosa) } 

(p x =d/dx[— ztana+k{-iztana(x+ztan 2 a/2)+(iMyztana/cosa)}] = 

2 

=— tana— iktana[x+ztan a-(Myz/cosa)]= 

2 

=-tana— ik[tana(x„+ztan a)— (y/ sinacosa)(sina/ cosa)]= 

=tana — ik[tana(x 0 + ztan 2 a)+(y/ cos 2 a)]= 

Therefore, the pressure distribution is given by the following relationship: 

2 

cp=2tana{ 1— ik[tana(x 0 + x (tan a-l)+2ycota} 

Using the results for zones I,I*,II JI* it is possible to calculate the pressure jump across 
the bottom airfoil (wing store airfoil) from the following relationship: 

Acp(x,0)=cp(x,0 + )- cp(x,CT) 

Two integration intervals, shown in Fig. 4, can be distinguished in the lower airfoil taking into 
consideration that 0.5<A<1: 

1) 0<x<A+B 

2) A+B<x<2A 

In the first integration interval: 

Acp I (x,0)=cp(x,0 + )— cp(x,0~)= 

2 

=2tana{ 1— ik[tana(x 0 +x (tan a-l)+2vcota} 

+2tana{ 1— ik[tana(xo+x (tan 2 a-l)+2ycota} and y=0 

2 

Acp,(x,0)=4tana{ 1— ik(X(,+ x(tan a— 1)} 



In the second integration interval: 
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+2tana[ 1 - 2e^+ik {2e , ^[x 0 +(x— B) (tan 2 a— 1 ) +2A]-(x 0 + x (tan 2 a— 1 ) +2yco ta) } ] 

= + 2tana[ 1 — 2e'^+ik {2e'^[xo+ (x— B) (tan 2 a— 1) +2A]— 2(x 0 +x(tan 2 a— 1 )) } ] since y=0 
i5 

But e =cos5+isin8 and the last expression can be written: 

=+2tana[l— 2(cos8+isin8)+ik{2 (cos8+isin8)[xo+(x— B)(tan 2 a— 1)+2A]— 2(x 0 +x(tan 2 a— 1))}] 

This expression is a complex equation with real and imaginary part In order to calculate the 
pitch-damping coefficient only the imaginary part of Acp is needed. Therefore one has: 

ImAcp I (x,0)=— 4tana(x 0 +x(tan 2 a— 1)) 
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I m A c p n (x,0)=+2tana[— 2sin8/k+2cos5[xo+(x— B)(tan 2 a— 1)+2A]— 2(xo+x(tan 2 a— 1))] 

The lift and moment coefficient can be calculated from the following relationship: 

Cl=J *Acpdx 
Cm=oJ 1 Acp(x-x 0 )dx 
By definition: 
c m=9[ Cm 6 +ik Cm©'] 

Qv(9' — C M 0'u = 

(l/ik)J A+B ImAcp I (x,0) (x-x^dx+A+sl 1 ImAcp n (x,0)(x-x Q )dx= 

J A+B- 4tana(x 0 +x(tan 2 a_l)(x_x Q )dx+A+af2tana[-2sin5/k+2cos5[xo+(x— B)(tan 2 a-1)+2A]— 2( 
Xo+x(tan 2 a-l))] (x-x^dx 

The last two integrals will be calculated separately. The first one is: 

2tana[-2xo 2 (A+B)+ x 0 (A+B)(A+B)-2(tan 2 a-l)[x 0 (A+B)72-(A+B) 3 /3)] (TV-1 2) 

The second one can be calculated using Maple. The results are a complicate expression. 
Two calculations have been included for the convenience of the reader since in both integrals 
there is a common factor 2tana. The result is: 

2tana{2x 0 +tan 2 a-B 3 A 2 B— cos8-tan 2 a*B 3 /3+2tan 2 a/3-sin8-A 2 /k— sin8-B 2 /k+sin8/k+2cos8/3 
+2xosin8-A/k— 2x 0 2 cos8-B-2x 0 2 cos8*A+4x 0 cos8-A+2x 0 cos8-B+2x 0 sin8-B/k+2x 0 2 A+2x 0 2 B 
+2cos8-tan 2 aA 3 /3+2cos8-AB 2 -4x 0 AB-2sin8-AB/k-2/3+4cos8-A 3 / 3+2B 3 /3+2AB 2 +2BA 2 
— 2tan 2 a-AB 2 -2tan 2 a-AB 2 -2cos8-tan 2 a/3— 2x 0 cos8-x 0 tan 2 a-cos8-B-2cos8-A-2x 0 sin8/k+2x 0 2 
sin8/k+2x 0 2 cos8-2x 0 2 +x 0 tan 2 a*A 2 +x 0 cos8-tan 2 a+x 0 tan 2 a*B 2 +cos8-Btan 2 a— 2xoCos8-Btan 2 a+ 
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3cosSA 2 B— 2tan 2 a*A 3 / 3— 2tan 2 oc*B 3 / 3+cos8-B 3 / 3+2A 3 /3— 2x 0 B 2 +2x 0 tan 2 oc-AB 
— 2 xqA 2 — 2x 0 cos5*BA+x 0 cos6*tan 2 a*B 2 — XqCosS-A 2 — XoCOsS-tan^-A 2 } (IV-13) 

Therefore: 

C M0 '=2tana[— 2 xq 2 (A+B)+ Xo(A+B)(A+B)-2(tan 2 a-l)[x 0 (A+B) 2 /2-(A+B) 3 /3)] + 
2tana{2x 0 +tan 2 a-B 3 A 2 B-cos5*tan 2 a-B 3 /3+2tan 2 a/3— sin5*A 2 /k— sin5*B 2 /k+sin8/k+2cos5/3 
+2x 0 sin5*A/k— 2x 0 2 cos5*B— 2x 0 2 cos5*A+4x 0 cos5*A+2x 0 cos8*B+2x 0 sin5-B/k+2x 0 2 A+2x 0 2 B 
+2cos8*tan 2 aA 3 /3+2cos8A.B 2 — - 4x^6— 2sin5*AB/k— 2/3+4cosS*A 3 /3+2B 3 /3+2AB 2 +2BA 2 
—2tan 2 aA.B 2 —2tan 2 a*AB 2 —2cos8*tan 2 a/3— 2x^058— x 0 tan 2 a-cos8-B—2cos8*A—2x 0 sin8/k+2x 0 2 
sin8/k+2x 0 2 cos8— 2x 0 2 H-x 0 tan 2 a*A 2 +x 0 cos8*tan 2 a+x 0 tan 2 a*B 2 +cos8*Btan 2 a— 2x 0 cos8-Btan 2 a+ 
3cos8*A 2 B— 2tan 2 a-A 3 / 3— 2tan 2 a*B 3 / 3+cos8*B 3 / 3+2A 3 / 3—2x 0 B 2 +2x 0 tan 2 a*AB 
— 2 xqA 2 — 2x 0 cosS-BA+x 0 cos8*tan 2 aB 2 — XqCosS-A 2 — x 0 cos8*tan 2 a*A 2 } (IV-14a) 

This is the major result of our investigation. It represents the pitch -damping coefficient of the 
lower airfoil in the presence of the upper airfoil. Note that both airfoils are oscillating in pitch. 
Also the two airfoil oscillation may lag by a phase angle 8. 

If B=0 and 8=180° the above result for the pitch damping coefficient is the same as 
that of an oscillating airfoil mounted close to a stationary airfoil or wall at a distance d/2(M. F. 
Platzer, H. G. Chalkey). As will be shown later the same result can be derived using the 
elementary theory. The pitch-damping coefficient for that case is given by the relation: 

CM6'- _ 2tana{x 0 2 (4-2 A)— Xq( 3A 2 -4A+2+(A 2 — 2)(tan 2 a-l))+2A 3 — 2A— 2(2— A 3 )/3 (tan 2 a-l)] 

(IV-14b) 



Where A— 2hcota, and h is the airfoil distance. 
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Fig. 5 



For the special case of zero stagger B=0 and out-of-phase oscillation d=180 the 
general result equations (IV-12) and (TV-13) simplifies to (TV-14). 
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This agrees with the result previously obtained by Platzer and Chalkey (1972) if 
A=d.cota in replaced by A=2hcota in replaced by A=2hcota used by Platzer and Chalkey. 
Two unstaggered airfoils oscillating out-of-phase are equivalent to a single mounted at a 
distance h=d/2 from a solid wall of a non-moving airfoil. 

It can be seen from Fig. 5 that the interference is highly destabilizing with increasing 
amounts of interference. 

Note that A=0.8 implies a larger amount of interference than A=0.9, for example. At A=0.9 
instability is encountered from low supersonic flow to a maximum of M=3, which at A=0.8 
instability occurs up to much larger supersonic Mach numbers. 




Another interesting case occurs for in-phase oscillation 8=0 of two unstaggered airfoils. The 
instability boundaries are shown in Fig. 6. Again, it is seen that increasing amount of 
interference leads to increasing regions of instability. 
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Fig. 6 shows curves of zero pitch damping coefficient (with torsional stability 
boundaries) for an airfoil mounted close to another airfoil when B=0 but oscillating in phase. 

Following the same procedure a special case of the above will be calculated. In that 
case the upper airfoil is assumed to be stationary while the lower airfoil (store airfoil) oscillates. 
The results will be compared with previous analyses (M.F. Platzer, 1971). 

Starting again with Sauers solution for the slowly oscillating airfoil in unbounded 
supersonic flow it is shown how this solution can be extended to consider the upper airfoil 
(wing) interference. To further simplify the comparison it is assumed that: 

4m 1 - 1 — < 1 

c 

Therefore, the moving airfoil is mounted at half distance from the stationary upper airfoil in 
compare to the previous analysis. The solution is given by the following set of equations: 

X ( x >y)=g(z) 

T (x,y) =h (z)+iMyg(z)/ cosa where z=x-ycota 
T— h(z)+iMy g(z)/cosa where z=x+ycota 

h(z) and g(z) are arbitrary functions for positive arguments of z and zero for z<0, it is obvious 
that the expressions of the above functions should be consistent with the last assumption. 

Using the two solutions for left and right running Mach waves in the supersonic flow 
field of slowly oscillating airfoils pressure coefficients can be obtained. In order to do so the 
flow field between the two airfoils should be divided into several zones. The number of zones 
depends upon A. It is obvious that for A>1, there is one zone along the airfoil, for 0.5<A<1 
there are two. So there are three regions. In the first zone I there is no interference from the 
airfoil wing (upper wing), which means that the lower airfoil (zone I only) does not sense the 
upper airfoiL In that case, the solution has the form: 

(PfcYrk) = g(z) +k{h(z)— iMyg(z)/ cosa} 
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The boundary conditions that must be satisfied at y=0 are: 

Xy =-cotcc-g'(z)=-l 

v Fy=-cota-h'(z)-iMyg(z)/ cosa=— i(x -Xq ) 

From the last two equations it is obvious that: 
g(z)=ztana 

h(z)=-iztana(x 0 +ztancx/2) 

Assuming a slowly oscillating airfoil, the solution to the PDE has the form: 
=g(z)+k[h(z)-iMyg(z)/cosa] 

When z = x— ycota 

In order to calculate the pressure distribution on the upper surface of the lower airfoil 
in the zone I the following relation should be used: 

cp =-2[ik(p +cp x ] 

Substituting the known functions cp and (p x 
cp=ztana +k{— iztana(x 0 +ztan 2 a/2)— ^Myztana/cosa)} 

(p r =d/ dx[ztana +k{-iztana(x 0 +ztan 2 a/2) — (iMyztana/cosa)}]= 

= tana-iktana[x 0 +ztan 2 a-(Myz/cosa)]= 

= tana-ik[tana(x 0 +ztan 2 a)+(y/sina-cosa)(sina/cosa) ]= 

=tana-ik[tana(x 0 +ztan 2 a)+ (y/ cos 2 a)] = 



45 



=tana-ik[tana(x 0 +(x-ycota)tan 2 a+(y/cos 2 a)] = 

=tana-ik[tana(x 0 + x tan 2 a-ytana+ (y/ cos 2 a)] = 

=tana-ik[tana(x 0 + x tan 2 a-ycota] 

Therefore, the pressure distribution is given by the following relationship: 

cp=-2{tana-ik[tana(x 0 +xtan 2 a-ycota]+ik[ztana+k{-i2tana(x 0 +ztan 2 a/2)- iMyztana/ 
cosa}]} 

Taking into consideration that higher order terms of k are neglected the pressure 
distribution is given by the following relationship: 

cp=2tana{— l+ik[xo+ x (tan 2 a-l)+2ycota]} 

So at y=0 

cp=2tana{— l+ik[x 0 + x (tan 2 a-l)]} 

In the zone II: 

9(>sy>k) =g(z) +k[h(z)-iM(yg(z)/cosa]+g 1 (zi)+k[h 1 (zi)+iM(yg 1 (z 1 )/cosa] 
where: 

zi=x-A+ycota 

In order to calculate the functions g(z),h(z) the boundary conditions will be applied 
The boundary conditions that must be satisfied at y=h are: 

X y =-cota-g'( z)+cota-g/( zj) = 0 
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t f 

'P y =-cota-h'(z)+cota-h 1 '(z 1 )-iM[g(z)-g 1 (z 1 )]/cosa+iMAcota[g (z)+g! (zj) ]/2cosa= 0 

After algebraic manipulation: 

cp— — 2tanoc {—1 4-ikfxo+x (tan 2 a— 1)]+2A— 2ycota} 

For the region III: 

<p(x,y,k) =g(z)+k(h(z)-iM(yg(z)/cosa]+g 1 (zi)+k[hi(z 1 )+ iM(yg 1 (z 1 /cosa]+ 

g 2 (z 2 )+k[h 2 (z 2 )-iMyg 2 (z 2 )/ cosa] 

where: 

z 2 =x-A-ycota 

The boundary conditions that must be satisfied at y— 0 are: 

Xy =- cota-g r ( z) +cota-g/( zj) -cota-g 2 '( z 2 )= -1 

'Fy=-cota-h'(z)+cota*h 1 , (zi)-cota-h 2 '(z 2 )-iM[g(z)-g 1 (z 1 )+g 2 (z 2 )]/cosa = - i(x -xo) 

After several manipulations: 

cp=— 2tana{3+ik[— 3xo+ (3x— 2A) (tan 2 a— l)]-2A/cos z a} 

Using the results from the three zones, the pressure jump across the lower blade can be 
calculated from the following relationship 

Acp (x,0) =cp (XjO 4 )— cp(x,0”) 

In order to calculate the cp(x,(T) previous results will be used (two airfoil oscillating case) for 
the lower surface of the store airfoil(store): 
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cp— 2tana { 1 +ik[x 0 +x(tan 2 a— 1)]} 



So Acp (x,0) — cp (x,0 + )— cp(x,(T) 

Two integration intervals can be distinguished in the lower airfoil taking into consideration that 
0.5<A<1: 

1) 0<x<A 

2) A<x<2A 

In the first integration interval is: 

Acp,(x,0)=cp(x 5 0 + )- cp(x,Cf)= 

=-2tana{ 1— ik[tana(x 0 + x(tan 2 a— 1)} = 

=— 2tana{ 1— ik[tana(x 0 + x (tan 2 a— 1)} and y=0 
Acpj(x,0)=4tana{l— ik(xo+ x (tan 2 a-l)} 

The latest expression is a complex number with real and imagin ary part. In order to calculate 
the pitch-damping coefficient only the imaginary part of Acp is needed. Therefore: 

ImAcpj(x,0)=4tana{— Xo+ x (tan 2 a-l)} 

In the second integration interval: 

ImA cp n (x,0) =+2iktana[— 4 xq+ (4x— 2A) (1 -tan 2 a)— 2A/ cos 2 a] 

The lift and moment coefficient can be calculated from: 

Cl=J 'Acpdx 

1 Acp(x-x 0 )dx 
By definition: 
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Qvt-Qol 0*9 +ik C^'] 



0*9' 0*9'l+ 0.l9'n 

For the purpose of the analysis, only the pitch-damping coefficient is required. Therefore, in 
order to facilitate the calculations only the necessary imaginary part of the will be 
calculated. 

0*9'= Q*e'i + Q*9'n = (1 / ik) J A ImAcp^XjO) (x-x Q )dx+ J 1 ImAcp n (x,0) (x-x^dx= 

J A — 4tana (xo+x(tan 2 a- 1) (x-x^cIx+aI 1 2iktana[-4xo+ (4x-2A) (1 -tan 2 a)-2 A/cos 2 a] (x-x^dx= 

=-2tana{x 0 2 (4— 2A)-x 0 (3A 2 -4A+2+(A 2 -2)(tan 2 a-l))+2A 3 -2A-2(2-A 3 )/3 (tan 2 a-l)] 

The above result is in agreement with the previous results derived by M. F. Platzer and 
H.G. Chalkey but the last one has been derived using the elementary theory. 

Therefore the elementary theory that has been developed predicts the above result as a 
special case for B=0 and 5=180° and represents for that case the pitch damping coefficient of 
an oscillating airfoil mounted close to a stationary airfoil or wall at a distance h/2. It is also the 
same result of the analysis of supersonic flow past a slowly two-dimensional airfoil in a wind 
tunnel with porous wall as presented by M.F. Platzer, 1971. 
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The second general case that will be examined in the thesis is shown in fig. 2. In that 
case, the Mach cone is wide enough to include the wing airfoil 




Using Sauer’s solution as previously and dividing the flowfield in zones as depicted in 
Fig.7 the perturbation potential (p in zone I is given: 

<P(x,y>k) =g(z)+k[h(z)+ iMyg(z)/cosa] 

Where: 

z=x-A-B +ycota at y=h 

The boundary conditions that must be satisfied at y=h are: 



Xy =cota-gf ( z) = -e 1 ^ 
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*8 

= cota-h'(z)+iMg(z)/cosa+iMy g'(z)/cosa =-ie l (x— xq-B) : 

After integration and algebraic manipulations: 

g(z) =-ztana-e 1 ^ 

c 

h(z) =iztana-e‘ (xo+ztan 2 a/2+Mh/ cosa) 

Where again the condition has been imposed that the functions g(z) and h(z) are zero for 
negative and zero arguments. 

*§ ‘S 

<t> r =-tana-e’ +iktana-e lO {(x 0 +z 1 tan 2 a+Mh/ cosa)-My/ cosa} 

= -tana-e 1 +iktana-e* {(xo+(x-A— B+ycota)tan 2 a+Mh/cosa)— My/cosa} 

*S ’8 

=tana-e* — iktana-e 1 {(xo+(x— B)tan 2 a— A tan 2 a+ycota+h/sina-cosa) 
y/ sina-cosa} = 

*S *8 

=-tana-e* +iktana-e‘ {(xo+(x— B)tan 2 a+A— ycota} 

Because h=A/ cota , tana-1 /cosa-sina=cota 

This gives the following result for the pressure distributions on the lower surface of the upper 
airfoil (wing airfoil) in zone I: 

cp =-2fikcp +<p x ] 

cp=2tana-e‘^— 2iktana-e‘^ { (xq+ (x-B) (tan 2 a— l)+2A-2ycota} — 
cp= 2tana-e‘^ { 1 — ik { (xo+(x— B) (tan 2 a-l)+2A— 2ycota} 
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cp= 2tana-e i5 {l-ik{(xo+(x-B)(tan 2 a-l)} 

In order to calculate the pressure on the upper surface of the wing airfoil the same approach 
will be followed and the result is: 

cp= 2tana*e 1 ^{-l+ik{(x 0 +(x-B)(tan 2 a— 1)} 

To calculate the Acp of the upper airfoil (wing airfoil) in the interval of zone I we have to 
calculate: 

Acp (x J h) = cp (xh*)-cp (xdO 
Therefore: 

Acp= 4tana-e’^ { 1 -ik { (xq+ (x-B) (tan 2 a-l)} 

Following the same arguments as previously for the zone II: 
zi=x-A— ycota -B 

gi(zi) =z 1 tana-(l-e‘ 5 ) 

The second boundary condition gives: 

-cota-h 1 , (z)+cota-h'(z)— iM[— g(z)+g 1 (z 1 )]/cosa=— i(x -xo— B-A) 

.C 

-cota-hi'(z)+icota*tanae lO (x 0 +ztan 2 a+Mh/cosa)-iMzitana/ cosa+ i(x - xo-B-A)=0 
Using the following trigonometric identity: 
ztan 2 a-z/cos 2 a=-z we have: 
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hi(Zi) =izitana{(xo+ z 1 tan 2 a/2)(-l+ e'^)+Mh e^/cosa}. 



In order to calculate the pressure distribution on the upper surface of the lower airfoil in the 
zone II the following relation should be used: 

cp =-2[ik(p +(pj 

Substituting the known functions (p and (p r 

(p=Zitana-(l-e 1 ^)— ztana*e , ^+k{iz 1 tana{x 0 +Zitan 2 (X/2)(-l +e l ^)+Mhe‘^/ cosa} +iztana-e’^(x 0 + 
ztan 2 a/2 +Mh/ cosa) } 

<p i =d/dx[z 1 tana- (l-e 1 ^)-ztana-e‘^+k{iz 1 tana{x 0 + zitan 2 a/2)(-l+e i ^)+Mhe 1 ^/cosa}+ iztana 
e‘^ (xq+ ztan 2 cx/2 +Mh/ cosa).}]= 

=tana-(l-e‘^)-tana*e , ^+k{itana{(xo+z 1 tan 2 a)(-l+e , ^)+Mhe i ^/cosa}+ itana-e’^(x 0 + ztan 2 a 
+Mh/ cosa)} 

= tana*(l— 2e l ^)+kitana[{(xo+ z ; tan 2 a)(— l+e'^+Mhe 1 ^/ cosa} +e‘^ (xo+ ztan 2 a+Mh/cosa)] 
T akin g into consideration that y=0 we have: 

=tana{ (1— 2e^)+ki[(xo+ Z!tan 2 a)(— 1 +2e‘^)+2Mhe'^/ cosa] } 

cp =— 2[ikcp +tp i ]=— 2[tana-(l— 2e‘^)+kitana{(x 0 + Zitan 2 a— Z[l)(— l+2e‘^)+ 2Mhe‘^/cosa}] 

The last relation gives the cp on the upper surface of the lower airfoil (store) in the zone ILIn 
the zone III a similar approach will be used to calculate the h 2 (z;>) and g 2 (z 2 ). The variable z 2 , 
taking into consideration the fig. 2, is given by the following relationship: 
z 2 = x— 2A+ycota — B 
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The first boundary condition that must be satisfied at y— h is: 



=z 2 tana-(l-e‘ ; 

The second boundary condition gives: 

-cota-h 1 '(z 1 )+cota-h'(z)+cota*h 2 '(z 2 )+iM[g(z)-g 1 (z 1 )+g 2 (z 2 )]/cosa+tMhcota[g'(z)+gi'(zi) 
+g 2 ' (Z 2 )] / cosa=-i(x-xo-B-2A) e‘^ 

cota-h 2 '(z 2 )+icota-tana(x 0 +z 2 tan 2 a)-e l ^Mz 2 tana/cosa-f-iMh(2— 3e^)/cosa4i(x— xo-B-2A)e^ 
=0 



cota-h 2 '(z 2 )=-i(x 0 +z 2 tan 2 a)+e I ^z 2 /cos 2 a— iMh(2— 3e^)/cosa— i(z 2 — Xo) e^ 

h 2 (z 2 )=— iz 2 tana (x 0 +z 2 tan 2 a/2)(l-e 1 ^)- iMh tana Z 2 (2— 3e*^)/cosa 
In a similar manner, as previously, the cp can be calculated: 

cp=z 1 tana-(l-e‘^)-ztana-e 1 ^+z 2 tana-(l-e l ^)+k{izitana{xo+z 1 tan 2 a/2)(-l+e' 5 )+Mhe' 5 /cosa} 

+iztana-e‘^ (xq+ ztan 2 a/2+ Mh/ cosa) +itana(x 0 + Z2tan 2 a/2)( 1 -e‘^)-iMhtanaz 2 (2-3e'^)/ cosa 
“ iM y[g( z )-gi( z i) + g 2 (z 2 )]/ cosa} 



cp x =tana-(2— 3e l ^)+kitana[(xo+z 1 tan“a)(— l+e‘^)+Mhe 1 ^/ cosa} +e i ^(x 0 +ztan 2 a+Mh/ cosa)+ 
(x^ z 2 tan 2 a)(-e l ^+l)-Mh(2-4e‘^)/ cosa] 

Cp I =tana < (2— 3e , ^)+kitana[(x 0 +Zitan“a) e 1 ^— Mh(2— 6e 1 ^)/ cosa] 

In order to calculate the pressure distribution on the upper surface of the lower airfoil in the 
zone II the following relation should be used: 
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cp=-2[ik(p +cp x ] 



Substituting the known functions cp and (p x 

cp=-2iktana-[(2-3e , ^)+z 2 (2-3e‘^)+ ( xq+ z 2 tan 2 a)e‘^-Mh (2-6e 1 ^)/ cosa] 



=-2iktana- [(2-3e , ^)+2z 2 (l-e 1 ^) + (xq+ (z 2 -l)tan 2 a)e^-Mh(2-6e‘^)/cosa] 
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V. METHOD OF CHARACTERISTICS 



In the method of characteristics, it is desirable to define the coordinate system in such 
a way that all-possible discontinuities could occur across. Along this coordinate system, the 
equations of motion of the flow field can then be treated as ordinary differential equations that 
are solvable by classical or numerical techniques. 

To obtain the equations in this coordinate system (the characteristic directions) in the 
(x,y) plane, the equations of motion are written in terms of the arbitrary intersecting 
coordinates: 

5 = 4(*>y) (v- 1 ) 

And 

*1 = V(x,y) (V- 2 ) 



If the first derivatives of the dependent variables, u,v,p with respect to £ are made 
indeterminate across lines of T) = constant, and the first derivatives of the dependent variables 
with respect to T| are made indeterminate across lines of %— constant, then any possible 
discontinuities in the first derivatives will occur across these lines. These lines are then the 
characteristics and their equations are obtained in the evaluation of the indeterminades. 
Consider first, two-dimensional, steady flow. The governing equations of motion are the 
continuity equation: 



d(P-v ) _ 0 

dx dy 



(V-3) 



The Euler equation becomes: 
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- du du 

u Hv — 

dx dy 



+±^.=0 

p dx 



— du dv 
u — + v — 
dx dy 



+ !^ = o 
p fy 



In addition, the energy equation becomes: 



- 3$ ds 
u — + v — 
dx dy 



= 0 



(V- 4) 



(V- 5) 



(V-6) 



Along a streamline: 



-dp dp -2-dp dp 1 

u — + v— -c [u— +v— 1 
dx dy dx dy 



(V-7) 



Substituting the above into the continuity equation one obtains: 



du dv 
dx dy 



-2 

C 



-dp dp 
u— + v— 
dx dy 



(V-8) 



In terms of the new coordinates: 



P£r«* + Pty v s +4-k +V ^K = ~P r h U n ~ P r ly V n “4*^ + v7 7yK 



In similar fashion the Euler equations become: 
/fe + V £y + &P{ = -pfo* + V7 7y K ~ Prfl* 
pfe + v£ y \> f + Z x p 4 = -pfyn x + V7] y ^ - p tJ rj x 



(V-9) 



(V-10) 

(V-ll) 
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The last three equations form a system of three equations Solving for 

by Cramer’s rule and taking into consideration that is indeterminate across T| both 
determinants in numerator and denominator should vanish. Therefore, from the denominator 
determinant we have: 

[< +v#,H(w 2 -~c)$ x 2 + 2v<^+(v 2 -c 2 )i/] = 0 (V-12) 



The solution to the above equation gives the equations of all three characteristics in the 
physical (x,y) plane: 



dx 



v 

— = tan^ 

u 



(V-13) 



dy 

— = tan (£±a) 

dx 



(V- 14) 



Where C, is the angle the streamlines makes with x-axis and a is the Mach angle. The first 
equation describes a streamline while the last one describes left and right running Mach lines. 

The compatibility equation for u,v,p along £ and ^characteristics can be obtained by setting 
the numerator of p^ equal to zero: 



du dv _ 1 dp 

v— — u — + — cot or— = 0 
drj dr] p dr] 



(V-15) 



The above relation must be satisfied along characteristics (£,T|). 

According to Teipel (1962) the unsteady flow over a flat plate can be treated in a similar 
fashion. In order to facilitate the calculations the continuity and Euler equations can be 
rewritten using the local sonic velocity: 

-2 p 

c = y— (V- 16) 

P 
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So canceling higher order teems the continuity equation becomes: 



2 dc 2 dc du dv 

■ + u„ — + c— + c— = 0 



y-\ dt y - 1 ” fix “ fix “fix 



(V- 17) 



In a s imilar , fashion the Euler equations become: 



du 2 dc du . 

— — H c x +W = 0 

dt y - 1 fix fix 



(V-18) 



fiv 2 5c dv n 

1- c m 1 - w_ — = 0 

dt y- 1 fix fiy 



(V-19) 



A system of three differential equations has been formed for three unknowns u,v,c 

To simplify the calculations with the assumption of harmonic oscillation we introduce the 
amplitude functions: 



u{x, y y = 



u_ 



(V-20) 



V{x,y)e ta> = 



1 v 
a/m 2 -1 



(V-21) 



C(x,y)e = 



icni 2 1 C 



Va / 2 -1 A / 2 c „ 



(V-22) 



Therefore, the continuity and Euler equations become: 



— + Vm 2 -1 — + M 2 —+ ikM 2 C = 0 
fix dy dx 



(V- 23) 



fix fix 



(V-24) 
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dx y]M 2 - 1 dy 



(V- 25) 



With the above variables the compatibility relations for the above system of equations are: 



r dV\ 

< dx) 4 



dC 

K dxj 4 



+ ik 



V + — 1 — (m 2 C-u) 
m 2 -v ’ 



= 0 



(V-26) 



(—) -(*) + ik 

te) n V&J, 



V — \ — (m 2 c-u) 

M 2 - V ’ 



= 0 



(V-27) 



r dU} 

& J str 



r dO 

<dx)*r 



+ ikU = 0 



(V- 28) 



Taking into consideration the irrotationality: 



du dv _ 
dy dx 



(V-29) 



The non-dimensional system of equations can be written as: 



~ + ylM 2 -\~ + M 2 — + ikM 2 C^0 
dx dy dx 



(V-30) 



6U 8C ;rr n 

+ — + ikU = 0 

dx dx 



(V-31) 






dy 



dx 



(V-32) 



Transforming the last system of equations to the new coordinate system (^,ri): 



'dU} fdV ^ 

. d* ) 4 l dx ) 4 



+ ik—¥——(U - C) = 0 
M 2 - V 



(V-33) 
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(V-34) 



r dU_ 

K dxy 



+ 



r av} 

U *J, 



+i ^( y - c )=° 




+ ikU = 0 



(V- 35) 



The last system of equations can be solved using finite differences. The values of U,V 
and C on the upper and lower surfaces of the airfoils can be calculated using the results of the 
Elementary Theory (Chapter III). The following relations apply: 



U = <p x ,V = tan a<py, C = -\cp x + ik<p\ 



(V-36) 



Using the results of chapter III the following results apply to zone I: 

Uj = ik{x 0 + xiax?" a + ycoXa)\ (V- 37) 



{- 



M 



V, = tana - 1 + ik{x - x + y(cota + )) 

1 1 0 cosar 



(V- 38) 



Cj = tanaj- 1 + + x(tan^ a - 1) + 2>>cota)]j (V-39) 



The last three equations represent the values of the parameters on the upper surface of the 
lower airfoil in the interval of the zone I. The last equation gives the pressure distribution and 
it is half the value of the actual pressure distribution coefficient. Following a similar approach 
the following results apply in the zone II* for the pressure coefficient: 



C = 2(2 tan ae^ - tana + ik tana[x Q + x(tan^ a - 1)] - 
2ike*$ tana[x 0 + (tan^ a - l)(x - B) + 2d cot a]} 



(V-40) 
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VI. FLUTTER ANALYSIS 



The lift and moment are given from the following relations: 



L = 



2 



pu>\c„dx 



0 



(VI-1) 



M = \pu 2 \c p {x-x 0 )dx 

4 0 

And 



(VI-2) 



Cp~-2J*£ +U &\ 

U 2 l dt 8x j 



(V 1-3) 



While 




(VI-4) 







And defining: 



c m=0 0 [c me + ikc m0 ,] 



(V 1-5) 



(VI-6) 



The last relations have been used in the elementary theory as shown in chapter III, and 
are a result of the following analysis. To compute the non dimensional lift and moment acting 
on the store airfoil the non-dimensional pressure coefficient P as defined in problem 
formulation last equation will be utilized and consequendy the lift and moment can be written 
as follows: 
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(VI-7) 



L = J{p(*,0 + )-P(;c,0-)}lt 

0 

M = J {P(*,0 + ) - P(x,0- )\x - x 0 )dx (VI-8) 

o 



And: 

P(x,y) = yM 2 C{x,y) (VI-9) 

Where: 

P(x,y)e ia * ? (VI-10) 

Pa, 

Although in the thesis only pitching, oscillation is considered the more general case 
include plunging oscillation also. In that case, Garrick and Rubinow’s (1946) method of 
expressing lift and pitching moment is used: 

+;I 2 ) + <? 0 (I, +/£„)] (VI-11) 

z c 

M = ~ Poo c Z uJk 2 e iat [^-(M l +iM 2 ) + 0 0 (M 3 + iM 4 )] (VI-12) 

4 c 



Where: 

k = — andc = 2b 
«- 

For a single degree of freedom as in that case the flutter analysis is much easier, and the 
equation of motion is: 



I + (1 + ig)C e 0 = M e 
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(VI-13) 



The sum of the moments of inertia about the elastic axis is: 

M f = -I e G (VI-14) 

In addition, the elastic restoring moment is: 

M R =-(l + ig)C e 0 (VI-15) 

Moreover, the aerodynamic pitching moment about the elastic axis is: 

M e = ~ P JcuJk 2 e'^0 o (M 3 +iM A ) (VI-16) 

Since: 

0 = 0 o e ,al (VI-18) 

0 = -(d 1 0 o e u * (VI-19) 

So the equation of motion becomes: 

e'“\.-l e m 2 9„ + (1 + ig)CA = + 1 M.) (VI-20) 

Separating the equation in real and imaginary parts gives: 

- I e o) 2 +C e +^pJ 4 0) 2 M, = 0 (VI-21) 

gC 0 + j pJ 4 co 2 M A = 0 (VI-22) 

Using the results of Garrick and Rubinow (1946): 

QgX - pr 2 +M 3 =0 (VI-23) 
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M a + gQgX = 0 



(VI-24) 



In order to solve the single degree of freedom problem it is necessary to obtain values 
of M 4 as a function of k until the second equation is satisfied The value of M 3 for that k is 
then used in the first equation to determine if X is reasonable: 



X = \- 



Ml 



>0 



(VI-25) 



The non-dimensional flutter frequency is: 
(Op 1 

(Oq -J ~X 



(VI-26) 



Accordingly the speed is: 

U F = 1 

C(O e kyfx 



(VI-27) 



Aerodynamic instability (flutter) occurs when M 4 < 0 while the L, M coefficients are given by 
the following: 



2 1 1 

L = L 3 + iL A = — r { |C(x,0 + )dx - JC(x,0" )a6c} 
k o o 



(VI-28) 



L 

M = M 3 +iM 4 = -- d-{JC(x,0 + )(x-x 0 )afr- JC(x,0")(x-x 0 )ofr} 

k 0 0 



(VI-29) 



And defining the following quantity, the non dimensional pitching moment, the analysis can 
be greatly facilitate: 

C m =d 0 [c m9 +ikc me .]e‘ h (VI-30) 
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(VI-31) 



c 



md 




(VI-32) 
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VII. CONCLUSIONS 



In this thesis, Platzer’s (1973) approach to analyze die aeroelasdc stability of flat-plate 
cascade and of airfoils in wind tunnels was extended to the case of two airfoils in close 
proximity to each other. The airfoils could oscillate in pitch with an arbitrary phase angle 
relative to each other. The analysis is valid for supersonic flow and is based on the small 
perturbation theory. 

Two •cases needed to be distinguished, namely airfoils with supersonic leading edge 
locus versus subsonic leading edge locus. In the former case, the leading edge of the upper 
airfoils is located upstream of the Mach wedge generated by the lower airfoil. In the latter case, 
the upper airfoil is within the Mach wedge of the lower airfoil 

For the case of supersonic leading edge locus a general formula could be derived for the 
pitch damping of the lower (store) airfoil which gives the dependence on Mach number, pitch 
axis location, phase angle, and proximity of the two airfoils. For the special case of out-of- 
phase oscillation between the two airfoils the result agrees with that derived by Platzer and 
Chalkey (1972), equations (37). The results showed that interference has a strongly 
destabilizing effect In addition, for the same special case Platzer’s (1973) equation (4.2) could 
be verified. 

Furthermore, the analysis of the subsonic leading edge problem was begun and the 
pressure distribution on both sides of the upper foil was derived. In addition, a method of 
characteristics procedure was described which could be used to compare results. 
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